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Abstract 

The cosmological baryon asymmetry can be explained as remnant of heavy Ma- 
jorana neutrino decays in the early universe. We study this out-of-equilibrium 
process by means of Kadanoff-Baym equations which are solved in a pertur- 
bative expansion. To leading order the problem is reduced to solving a set of 
Boltzmann equations for distribution functions. 



The generation of the cosmological matter-antimatter asymmetry in an expanding 
universe requires baryon number violation, C and CP violation, and a deviation from 
thermal equilibrium The classic mechanism which realizes these conditions is 
the decay of weakly interacting massive particles in a thermal bath 0. Particularly 
successful is the leptogenesis scenario where the decaying particles are heavy Majorana 
neutrinos 0]. The resulting baryon asymmetry is entirely determined by neutrino 
properties. The observed order of magnitude can be naturally explained without any 
fine tuning of parameters and in accord with present experimental indications for 
neutrino masses @]. 

The generation of a baryon asymmetry is an out-of-equilibrium process which is 
generally treated by means of Boltzmann equations. A thorough discription of the 
basic ideas can be found in 0. Some subtleties have recently been discussed in 0. A 
shortcoming of this approach is that the Boltzmann equations are classical equations 
for the time evolution of phase space distribution functions. On the contrary, the 
involved collision terms are S'-matrix elements which involve quantum interferences 
of different amplitudes in a crucial manner. Clearly, a full quantum mechanical treat- 
ment is highly desirable. It is also required in order to justify the use of Boltzmann 
equations and to determine the size of corrections. 

All information about the time evolution of a system is contained in the time 
dependence of its Green functions which can be determined by means of Dyson- 
Schwinger equations. Originally these techniques were developed for non-relativistic 
many-body problems. More recently, they have also been applied to transport phe- 
nomena in nuclear matter the electroweak plasma |TD|,|TI| and the QCD plasma 



r2| . Alternatively, one may study the time evolution of density matrices |[13|Jl4 
In the following we shall investigate non-equilibrium Green functions which are rele- 
vant for leptogenesis. We shall construct a perturbative solution of the corresponding 
Kadanoff-Baym equations which, to leading order, turn out to be equivalent to a set of 
Boltzmann equations. Higher-order corrections can then be systematically evaluated. 

Consider now the standard model with three additional right-handed neutrinos 
whose interactions are described by the lagrangian, 

£ = Il0AV^ - ]^v'^rMpr + h.c. (1) 

Here II and denote lepton and Higgs doublets, respectively. We shall restrict our 
discussion to the case of hierarchical Majorana neutrino masses. Mi <^ M2, M3. The 
baryon asymmetry will then be determined by the CP violating decays of the lightest 
Majorana neutrino Ni = vm + v'f^^ = N, 

r(Ar^ = 1(1 + e)r, r{N^U) = ^{i-e)r . (2) 

Here F is the total decay width and the parameter e <^ 1 measures the amount of 
CP violation. The generation of the baryon asymmetry takes place at a temperature 
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T ~ Ml = M -C M2,M3. It is therefore convenient to describe the system by an 
effective lagrangian where the two heavier neutrinos have been integrated out, 

L = ~lLi^\*^N + N^X-aChi(t)-\MN^CN 

+ \mjlu<P C + \riULi4> C , (3) 
with ^ 

For leptogenesis one has to consider the phase space distributions for heavy neu- 
trinos (/at), leptons (/;), anti-leptons (/[), Higgs (/</,) and anti-Higgs bosons (Z^). 
The generation of the lepton asymmetry is a process close to equihbrium. Hence 
one can hnearize the Boltzmann equations in the deviations from the equihbrium 
distributions. Due to the interactions in @ one has Sfi = —Sfj = 5f^ = —5f^. The 
Boltzmann equation for the Majorana neutrino reads 

d 

gN-g^SfN{t,p) = (5) 
-9N^^fN{p) - ^ / d^u{p)6fM{t,p) {\M{N ^ + \M{N ^ . 
For the lepton doublets one obtains 

2gipMt,k) = 
1 

2k 

1 

~2k 
1 

~2k 
1 

~2k 
1 

~4k 
1 

~2k 



d^l,{k) e 6Mt,Pi) {\M{N ^ + \M{N ^ U)]' 

d<l>,-,{k) {Sfi{t,k)U{p,) + fi{k)6U{t,p,)) 

X {\M{14> N)\^ + IMiU N)\^) 
d<l>m{k) {6 flit, k)U{p,) + mSU{t,pi)) 
X {\Mil(t> ^ r0)p + \Mm ^ 

d^-mik) {5fi{t,p,)U{p2) + fi{pi)5U{t,p2)) 

X [\M{l<p 10)1' + \M{1^ ^ 
d^mik) {5 hit, k)fi{p,) + fi{k)6Mt,p,)) 

X {\M{ll 00) p + \M{ll^4>^)\ 
d<!>m{k)6Mt,p,)U{p2)(\M{<p<P^ll)\' + \M{^^^ll)n . (6) 



Here 



<;*■...»... W - . . . . . . (2.) V(p + p. + . . . - - . . .) , (7) 

/.(p) = exp (-PEtip)) , (8) 

denote phase space integrations and distribution functions, respectively. The temper- 
ature T = and A4{. . .) is the matrix element of the indicated process, = '2 
and gi = gj = Q are the number of 'internal' degrees of freedom for the Majorana 
neutrino and the lepton doublets for three generations, respectively. For simplicity 
we have assumed small number densities so that we can use Boltzmann distribution 
functions for bosons and fermions and also neglect distribution functions for particles 
in the final state. The effect of the Hubble expansion is included by introducing 
the 'covariant' derivative d/dt — > d/dt — Hpd/dp. Integration over momenta then 
yields the more familiar form of the Boltzmann equations for the number densities. 
Eq. (H) describes the decay of the heavy Majorana neutrinos. Note, that also the 
equilibrium distributions are time dependent since the temperature varies with 
time. For massless particles the distribution functions are constant with respect to 
the 'covariant' time derivative. The first term in eq. drives the generation of a 
lepton asymmetry; the remaining terms tend to wash out an existing asymmetry. 
Eqs. d^) and (H) determine 6/^ and 6fi as function of time. We have only kept the 
interactions given by the lagrangian (0). A complete discussion can be found in |T3]. 



Green functions near thermal equilibrium 



The time evolution of an arbitrary multi-particle lepton-Higgs system can be 
studied by means of the Green functions of lepton and Higgs fields. For the heavy 
Majorana neutrino one has 

iGal3{Xi,X2) =TT{pTNa{xi)Nf^{x2)) , (9) 

where T denotes the time ordering, p is the density matrix of the system, the trace 
extends over all states, and the time coordinates ti and t2 lie on an appropriately 
chosen contour C in the complex plane |T^. G{xi, X2) can be written as a sum of two 
parts, 

G{xi, X2) = e(ti - t2)G>ixuX2) + e{t2 - ti)G<{xuX2) , (10) 

where 



iG>{Xi,X2)af3 = T^T^ {pNa{Xi)N/3{x2)) , iG'^ {Xi, X2) af3 = {pNp{x2) N^iXij) . (11) 

The 'time ordering' in eq. (|T0]) is along the contour C. 

For a system in thermal equilibrium at a temperature T = 1/(3 the density matrix 
is p = exp {—jSH), where H is the Hamilton operator. In this case the Green function 
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only depends on the difference of coordinates and it is convenient to introduce the 
Fourier transform, 

G(p) = J d^xe'P'^Gix) . (12) 

The contour C can be chosen as a sum of two branches, C = Ci U C2, which he above 
and below the real axis. The time coordinates are real and associated with one of the 
two branches. Correspondingly, the Green function becomes a 2 x 2 matrix, 

= [ G-(p) G-(p) j ■ 
The off-diagonal terms are given by 

G''{p)=G<ip), G'\p) = G>{p). (14) 



The diagonal terms of the matrix (|T3D are the familiar causal and anti-causal Green 
functions. The functions G^{p) and G^{p) satisfy the KMS-condition, 

G<{p) = -e-^P^G>{p) , (15) 

and the free Green functions are explicitly given by 

zG>(p) = {Q{p,)-Q{p,)f^{E)-Q{-po)ffjmPN{p). (16) 
iG<{p) = {Q{-p,)-Q{p,)f^{E)-Q{-po)h{E))p^{p) , (17) 

with the spectral density 

p^{p) = 2n{i + M)G-H{p' - , (18) 

and the Fermi-Dirac distribution functions 
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f^[E) = ff,{E) = -^^, E=^M^ + pK (19) 

Since N{x) is a Majorana field one has = fj^, and in the spectral density ( |T^ ) the 
charge conjugation matrix G occurs. In the following we shall also need the retarded 
and advanced Green functions, 

G^{x) = ±0(±x°) {G>{x) - G<{x)) , (20) 

which can be written as sum of an on-shell and an off-shell contribution, 

N_ , l/^>/ N ^^ , 1 ^ f J, ,fG>{x,u;',p) -G<{x,uj',p) 



G^ip) = ±- (G>(p) - G<ip)) + —V / duj' ^ ' _ (21) 

Z Zm J LU — u 

The Green functions for the lepton doublets and for the Higgs doublet, 

iS{xi,X2)ai3S^ = Tr [pTr^{xi)kfs{x2)) , zA(xi,X2)5," = Tr (pTr(xi)0,*(x2)) , (22) 



have the same structure as G{xi,X2)- The corresponding equations for S{p) are 
obtained from eqs. (p!0D -(p^ by replacing the spectral density Pn{p) by 

piip) = 2nPLMp' - M^) , P^ = i^, (23) 

and the distribution functions fN^E) and fN^E) by 

f^(E,p,) = fjiE,-f,i) = -^^^^^, E=\p\, (24) 

where ni is the lepton chemical potential. For the Higgs field one has 

p^{p) = 2n6{p^-M^) , (25) 

UiE, /i^) = f^{E, -fi^) = ^^(^_^^) _ ^ , ^=1^1- (26) 

We are considering a process close to equilibrium. This suggests that the cor- 
responding deviations of the Green functions may be obtained from the equilibrium 
Green functions by a small change of the distribution functions, 

i5G{x,p) = -6fNix,p)pNip) , (27) 

iSS{x, k) = -e{ko)Sfi{x, k)pi{k) , iSA{x, q) = -t{qo)5f^{x, q)p<p{q) ■ (28) 
Here we have used that due to the interactions given in (H) 6fi = —Sfj = 6f^ = —5f^. 

Kadanoff-Baym equations 

The Green functions for the heavy neutrino and the leptons satisfy Dyson- 
Schwinger equations, 

Ciidi — M)G{xi,X2) = 5{xi — X2) + / (i'^X3S(xi, X3)G(a;3, X2) , (29) 

- ..-.,./..n,.,.,.,.,.,, ,30, 

Jc 

where E and 11 are the corresponding self energies and the time integration is carried 
out along the contour C. Eqs. (^) and ( PP] ) can be turned into matrix equations 
with real time integration in the usual manner. For the off-diagonal elements and 
one then obtains 

C{i^i-M)G>{xi,X2) = j (fx^(T.>{xi,x^)G-{x^,X2) 

+S+(xi,a;3)G>(a;3,X2)) , (31) 

r^.S>ix..X2) ^ /d%(n>(.„X3)^-(x3,..) 

+n+(xi,a;3)5>(x3,X2)) . (32) 
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Here and 11^ are off-diagonal matrix elements of S and 11, which are defined 
analogous to . The equations for and S*^ can be obtained from eqs. (|3l| ) and 
(|32|) by replacing the superscripts '>' by '<'. Equations of the type (0), (|32D have 
first been obtained by Kadanoff and Baym for non-relativistic many-body systems 

i- 

For processes where the overall time evolution is slow compared to relative motions 
the Kadanoff-Baym equations can be solved in a derivative expansion. One considers 
the Wigner transform for G(xi,X2), 

G{x,p) = J d'y G + |,a; - I) , (33) 

and S{x,p), A{x,p), respectively. For the Wigner transform of a convolution one has 
in general, 

d^ye'Py J d^X2 A{xi,X2)B{x2,X3) (34) 

i I d d \ 

= A{x,p)B{x,p) - - \—A{x,p)—B{x,p) - —B{x,p)—A{x,p)j + ... , 

where x = {xi + X3)/2 and y = {xi ~ X3)/2. 

Using the derivative expansion (Q) the Kadanoff-Baym equations (0) and (|3^ ) 
become local in the space-time coordinate x. Keeping to zeroth order only the on- 
shell part of retarded and advanced Green functions and self-energies, which are given 
by expressions analogous to eq. (|2l|) , one obtains the equations 

GC-^ + ^-M)G>ix,p) = gC-^ + ^-M)G<{x,p) 

= ^i^>{x,p)G<{x,p)-^<{x,p)G>{x,p)) , (35) 

C-<^ + ^)S>{x,k) = C-^ + mHx,k) 

= ^(U>{x,k)S<{x,k) -U<{x,k)S>{x,k)) . (36) 

Solutions of these equations yield the first terms for the non-equilibrium Green func- 
tions G^{x,p)...S^{x, k) in an expansion involving off-shell effects and space-time 
variations, which include 'memory effects'. As an example for the type of corrections 
we list the first derivative term on the right-hand side of eq. (|35|), 

Ae = -l(|-S>(x,p)|-G<(x,j9)-|-S>(x,p)|-G<(x,p) 
A \ox op op ox 

-^S<(x,p)|g>(x,j>) + |s<(x,j>)^G>(x,j>)) . (37) 
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Figure 1: One-loop self energies for the Majorana neutrino. 



Solutions of the Kadanoff-Baym equations can be studied once the self-energies 
E and 11 are known. For weak coupling these can be determined in perturbation 
theory. 

Self-energies for lepton fields 

The one-loop contributions to the self-energy of the Majorana neutrino are shown 
in fig. (1). For vanishing chemical potential fig. (la), for instance, yields the result 



- ^Sg)>(p) = 2(AtA)n / 7^7^(2^)'<^(^^ -^^1 - P2)CS> {p^)A> {p,) . (38) 



In the following we shall only consider the case M ^ T, where the heavy neutrinos 
are non-relativistic. In this case the number density is small, /iv(-E') <S 1, a-nd one 
finds for the difference of the self-energies (po > 0) , 



-i(E>(p)-E<(p)) = -2(AtA)n /rf$i2(p)Cj^i 

X ((1 - fiiPim + fM)) + fi{Pi)U{P2)) (39) 



Here F is the vacuum decay rate of the Majorana neutrino. 

The one- and two-loop contributions to the lepton self-energy are shown in 
fig. (2a)-(2d). In the following we only list the terms which are needed for the solution 
of the Kadanoff-Baym equations to leading order. 

Particularly interesting are the terms fig. (2b) which drive the generation of an 
asymmetry. After some algebra one finds {k^ > 0), 



Here we have only given the deviation from the equilibrium self-energy which is ob- 
tained by using for the Majorana neutrino propagator the deviation from the equi- 
librium propagator 6G oc 6fN- Furthermore, we have again considered the case of 




(40) 



t (6U^'^<{k) + 6U^'>{-k)) = |-Im(Atr/A*)nM J d<^r2{k)^2PLSfN{t,P2) ■ (41) 
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Figure 2: One- and two-loop self energies for the lepton doublet. 
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small densities. The combination of Yukawa coulings is precisely the one occuring in 
the CP asymmetry of the Majorana neutrino decay. Since the heavy neutrinos N2 
and have been integrated out, the contributions fig. (2b) involve both, self-energy 
and vertex corrections ||T7| , p!8| , |T9| which, up to a numerical factor, are identical in this 
limit. 

All graphs in fig. (2) contain washout processes. In the case of small densities one 
obtains from fig. (2a) (/cq > 0), 

- z (ng)>(A;) - n(",)<(A;)) = -2(AtA)n / d^,^{k)i>,PLUM) + /^(P2)) • (42) 

For T ^ M the dominant contribution to the washout processes is due to fig. (2c) 
where the Majorana neutrino propagator can be replaced by a local interaction. For 
small densities one finds (fco > 0), 

- ^ (n(f (A:) - n(f (fc)) = -6^^ / rf$i23(^) {MiiPi) + 2^2/0(Pi)) Pl . (43) 
The complete expressions for the self-energies will be given elsewhere. 



Kinetic equations 



Given the lepton self-energies we can now look for solutions of the Kadanoff-Baym 
equations ( pSj) and (^61). A straightforward calculation shows that the right-hand side 
of these equations vanishes for equilibrium Green functions and self-energies. Since 
baryogenesis is a process close to thermal equilibrium we can search for solutions 
which are linear in the deviations, 

5G{t,p) = G>{t,p)-G>^{p) = G<{t,p)-G<{p), (44) 
6S{t,p) = S>{t,p)-S>{p) = S<it,p)-S<{p) . (45) 

One then obtains for the perturbations 5G{t,p) and 6S{t,p), 

zG^'pG{t,p) = zC^'^G>^{p) + {J:>^{p)-^<^{p))6G{t,p), (46) 

t^'pS{t,k) = (n>(fc)-n<(fc))<55(t,fc) 

+5n>(t, k)S<{k) - 5n<(t, k)S^g{k) . (47) 

The Green functions depend on time explicitly, as well as implicitly through the time- 
dependence of the temperature. Once the 'covariant' time derivative is used, the later 
vanishes for equilibrium Green functions of massless fields. This is not the case for 
massive fields. Hence, the first term on the right-hand side of (^Bf ) drives the deviation 
from thermal equilibrium. 
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We can now insert the perturbative expressions for the self-energies into eqs. (E6| 



and check whether the ansatz (pTl), (|28|) for 5G and 5S yields a solution. After 



some algebra one finds that this is indeed the case provided the distribution functions 
5fN and 5fi satisfy the following ordinary differential equations, 

E-5Mt,p) = -E- f ^(p) -2i\^\)n I d^r2{p)5fN{t,p)p- Pi, (48) 

gik^Sfi{t,k) = ^luiiX'^riX^^M j d^i2{k)5fN{t,p^)k-p^ 
ot 87r J 



-2(AtA)n j d^^^{k) (6 flit, k)U{p{) + fiik)6Uit,p,)) k ■ p^ 
-6^^^ / ^*i23(fc)(2(5//(t,A:)/^(pi) + /Kfc)5^(t,Pi) 

+Sfl{t,P2)U{P3) + fl{P2)SU{t,P3))k-p2 
H^flit,k)fM + m5h{t,p^) 

+2SU{t,p2)U{p3))k-pi) . (49) 

Comparing these equations with the Boltzmann equations (^) and (^ one finds that 
the two sets of equations are identical to leading order in the coupling where matrix 
elements and CP asymmetry are given by 

\M{N{p) ^ l{p^)<f){p2))\' = 4(AtA)nP-Pi, (50) 

\Mm^iPi)-KP2mP3))\' = 24^^^fc-p2, (51) 



3 Im(A^r/A* 



-M . (52) 



1677 (AtA)ii 

We conclude that for non-relativistic heavy neutrinos a solution of the Boltzmann 
equations generates a solution of the full Kadanoff-Baym equations to leading order in 
the expansion described above. For relativistic heavy neutrinos the matrix structure 
of the equations is more complicated and the time evolution of the different poles of 
the Majorana neutrino propagator are described by different equations. 

Given a solution of the Kadanoff-Baym equations to leading order the various 
corrections can be systematically studied. Note, that the solutions of eqs. (^Sf ) and 
(|49|) are not of the form Sfi(t,p) = hi(t)fi{p). Hence, the usual assumption of kinetic 
equilibrium does not appear to be justified. The size of 'derivative terms', which cor- 
respond to memory effects, and off-shell corrections can be determined by inserting 
the leading order solution into the various correction terms described above. Partic- 
ularly interesting are relativistic corrections in the case that leptogenesis takes place 
at temperatures T ~ M. 

The analysis of the Kadanoff-Baym equations for leptogenesis can be used to 
obtain constraints on the parameters M, (A^A)ii and e, which provides a quantitative 
relation between the cosmological baryon asymmetry and neutrino properties. 
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